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We calculate exactly by field theoretical techniques the sec- 
ond topological moment (m 2 ) of entanglement of two closed 
polymers Pi and P2. This result is used to estimate approx- 
imately the mean square average of the linking number of a 
polymer Pi in solution with other polymers. 

1. Consider two closed polymers P\ and P 2 which 
statistically can be linked with each other any number 
of times m = 0, 1, 2, . . . . The situation is illustrated in 
Fig. [I] for m = 2. 

Pi 



<m 2 > 



/ d 3 x\d 3 x 2 f_™ dm m 2 G m (xi, x 2 ; Li, L 2 ) 
J d 3 x 1 d 3 x 2 f*™ dm G m (xi, x 2 ; Li, L 2 ) 
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FIG. 1. Closed polymers Pi, P2 with trajectories Ci, C2 re- 
spectively. 

An important physical quantity is the probability dis- 
tribution of the linking number m as a function of the 
lengths of Pi and P 2 . As a first step towards finding 
it we calculate, in this note, an exact expression for the 
second moment of the distribution, (m 2 ). 

An approximate result for this quantity was obtained 
before in Ref. [|l| on the basis of a a mean-field method, 
considering the density of bond vectors of P 2 as Gaussian 
random variables. Such methods are usually quite accu- 
rate when a large number of polymers is involved ^,|). 
As an unpleasant feature, however, they they introduce a 
dependence on the source of Gaussian noise, and modify 
the critical behavior of the system, whereas topological 
interactions are not expected to do that [|,|[. Our note 
goes therefore an important step beyond this approxi- 
mation. It treats the two-polymer problem exactly, and 
contains an application to the topological entanglement 
in diluted solutions. The relevance of the two-polymer 
systems to such systems was emphasized in Q . Focusing 
attention upon a particular molecule, Pi, one may imag- 
ine all others to form a single long effective molecule P 2 . 

2. Let G m (xi, x 2 ; Li, L 2 ) be the configurational prob- 
ability to find the polymer Pi of length Li with fixed 
coinciding end points at xi and the polymer P 2 of length 
L 2 with fixed coinciding end points at x 2 , topologically 
entangled with a Gaussian linking number m. 

The second moment (m 2 ) is defined by the ratio of 
integrals B 



The denominator in (|l|) plays the role of a partition 
function: 

/r+00 
d 3 x x d 3 x 2 I draG m (xi,x 2 ;Li,L 2 ). ( 2 ) 
J — 00 

Due to the translational invariance of the system, the 
probabilities depend only on the differences between the 
end point coordinates: 

Gm (xi,x 2 ;Li,L 2 ) = G m (xi - x 2 ;Li,L 2 ) . (3) 

Thus, after a shift of variables, the spatial double inte- 
grals in ([l]) can be rewritten as 

J d 3 xid 3 x 2 G m (^i-^2]L 1 ,L 2 ) = V J d 3 xG m (x; L 1: L 2 ), 

where V denotes the total volume of the system. 

3. The most efficient way of describing the statistical 
fluctuations of the polymers Pi and P 2 is by two com- 
plex polymer fields ^^(xi) and ^2 2 (x 2 ) with n\ and n 2 
replicas (a\ = 1, . . . , ni, a 2 = 1, . . . , n 2 ). At the end we 
shall take m, n 2 — > to ensure that these fields describe 
only one polymer each 0. 

For these fields we define an auxiliary probability 
G\ (xi , x 2 ; z) to find the polymer Pi with open ends at 
xi , x^ and the polymer P 2 with open ends at x 2 , xj> . The 
double vectors x\ = (xi,x^_) and x 2 = (x 2 ,x 2 ) collect 
initial and final endpoints of the two polymers Pi and 
P 2 . Here we follow the approach of Edwards [|| , in which 
one starts with open polymers with fixed ends. The case 
of closed polymers, where m becomes a true topological 
number and it is thus relevant in the present context, is 
recovered in the limit of coinciding extrema. We notice 
that in this way one introduces in the configurational 
probability an artificial dependence on the fixed points 
xi and x 2 . In physical situations, however, the fluctua- 
tions of the polymers are entirely free. For this reason 
we have averaged in (|l|) over all possible fixed points by 
means of the integrations in d 3 xid 3 x 2 . 

The auxiliary probability G\ (x\ , x 2 ; z) is given by a 
functional integral H 



Ga(xi,x 2 ;z)= lim / P(fields) 

ni,ra,2— >0 J 

x (xi w w 2 2 (x 2 w 2 a2 (A)e~ A , 



(4) 
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where D(fields) indicates the measure of functional in- 
tegration and ai,<22 are now fixed replica indices. A = 
Acs + -4poi is the action governing the fluctuations. It 
consists of a polymer action 



where 



*4o =Acs+ d\ 



(10) 
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(5) 



and a Chern- Simons action to describe the linking num- 
ber m 



From Eq. ( |To| ) it is clear that Ga=o (xi, X2; z) is the prod- 
uct of the configurational probabilites of two free poly- 
mers. In fact, the fields ^2,^2 are f ree > whereas the 
fields ^1,^1 are apparently not free because of the cou- 
plings with the Chern- Simons fields through the covari- 
ant derivative D 1 . This is, however, an illusion: Inte- 



^4 CS =in [ d 3 xe^ vp A^d v A p 2 , (6) grating out A^ in (|), we find the flatness condition: 

J E^ vp d,A l . = 0. On a flat soace with vanishing" boundarv 



In Eq. (H) we have omitted a gauge fixing term, which en- 
forces the Lorentz gauge. The effects of self-entanglement 
and of the so-called excluded- volume interactions are ig- 
nored. The Chern-Simons fields are coupled to the poly- 
mer fields by the covariant derivatives = V + 27^ A*, 
with the coupling constants 7^2 given by 71 = ft, 72 = A. 
The square masses of the polymer fields are given by 
mf = 2Mzi, where M = 3/a, with a being the length of 
the polymer links, and Z{ the chemical potentials of the 
polymers, measured in units of the temperature. The 
chemical potentials are conjugate variables to the length 
parameters L\ and L2, respectively. The symbols col- 
lect the replicas ip^ of the two polymer fields. Let us note 
that in the topological Landau-Ginzburg model (||) the 
Chern-Simons fields do not change the critical behavior 
of the system, as expected. 

The parameter A is conjugate to the linking number 
m. We can therefore calculate the desired probability 
G m (xi, X2; Li, L2) from the auxiliary one Ga(xi,X2;z) 
by Laplace integrals over z = (zi,z\) and an inverse 
Fourier transformation over A. 

4. Let us use the polymer field theory to calculate the 
partition function (||) . It is given by the integral over the 
auxiliary probabilities 



c+00 



1 7 -< 1 , Mdz\ Mdzo ~ j j 

Z= \ d 6 x 1 d 6 x 2 lim / £^1+^2 

2m 2m 



dm / ^Ae- imA G A (xi,x 2 ;z). (7) 

-co J — CO 

The integration over m is trivial and gives 27r<5(A), en- 
forcing A = 0, so that 



/ 



Z = d 6 x 1 d 6 x 2 lim 



c+ico 



MdziMdz2 Li+Z2L2 



1 J c— too 



2m 2m 



x G\=o (xi,x 2 ;z) . 



(8) 



To calculate G\=o (xi , X2 ; z) we observe that the action 
A is quadratic in A. A trivial calculation gives 



G A =o (xi,x 2 ;z) = J £>(fields)e 



-Ao 



(9) 



e^ vp d v A l = 0. On a flat space with vanishing boundary 
conditions at infinity this implies A^ = 0. As a conse- 
quence, the functional integral (||) factorizes 

G\ =0 (xi,x 2 ;F) =G (xi -x / 1 ;zi)G (x 2 -x 2 ;z 2 ), (11) 

where Go(x^ — x^; Zi) are the free correlation functions of 
the polymer fields 



Go(*i 



A; Li) 



c+ico 



M 
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Thus we obtain for 



c— 100 
3/2 



Mdzj 
2m 

M(x i -x' i ) 2 /2L 



e ZiLi G (xi -xj;*0 

(12) 



the integral 



Z = 2tt d 3 xid 3 X2 lim Go(xi — x^; L\) Go(x 2 — x 2 ; L 2 ), 
J x i^ x i 



yielding the partition function 
2nM 3 V 2 



Z = 



(4tt) 3 



-3/2 



(13) 



It is important to realize that in Eq. (j?]) the lim- 
its of coinciding end points x.[ — > x^ and the inverse 
Laplace transformations do not commute unless a proper 
renormalization scheme is chosen to eliminate the diver- 
gences caused by the insertion of the composite operators 
|C(x)| 2 and |^(x)| 2 . 

5. Let us now turn to the numerator in Eq. (Eh. Ex- 
ploiting the identity m 2 e~ %rnX = —d 2 e~ %mX /d)?, and 
performing two partial integrations in A, the same tech- 
nique used above to evaluate the partition function Z 
yields 



N = K 



c+ico 

d 6 x\d 6 X2 lim 

2 ^0 ^C— ZOO 



Mdz x Mdz 2 
2m' 2m 



) ZlLi+Z2L2 



©(fields) exp(-A)|V'? 1 (xi)r^ 2 02 (x2)| z 
J d 3 xA 1 • **V*i^ + i J d 3 xA 2 1 |^i 



(/ 



d 3 xA 2 -^* 2 V^ 2 



X - j d 3 xA 2 2 \V 2 



(14) 
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where Ao has been defined in jTo|). In the above equa- 
tion we have taken the limits of coinciding endpoints in- 
side the Laplace integral over z\,z 2 . This will be justified 
later on the grounds that the potentially dangerous Feyn- 
man diagrams containing the insertions of operations like 
|^i| 2 vanish in the limit ni,n 2 — > 0. The functional in- 
tegral in Eq. ( |l4| ) can be calculated exactly by diagram- 
matic methods since only four diagrams shown in Fig. (||) 
contribute. 



FIG. 2. Four diagrams contributing in Eq. The lines 

indicate correlation functions of ^-fields. The crossed circles 
with label i denote the insertion of (~Xi)\ 2 . 

Only the first diagram in Fig. || is divergent from the 
loop integral formed by two correlation functions of the 
vector field. This infinity may be absorbed in the four- 
\I/ interaction accounting for the excluded volume effect 
which we do not consider at the moment. No diver- 
gence arises from the insertion of the composite fields 

IVf(x;)| 2 . 

6. In this section we evaluate the first term appearing 
in the right hand side of Eq. ( |l4| ) : 



— hm / 

n J c-l 



c+ioo 



J d 3 x\d 3 x 2 J d 3 x / 1 d 3 x / 2 



Md Zl Mdz2 eZlLl+Z2L2 
2ni 2ni 



(15) 



K 1 (xi)| 2 |V- 2 ° 2 (x 2 )| 2 (|*i| 2 Af) x/ (|* 2 | 2 Ai) x , ). 



There is an ultraviolet-divergence which must be regu- 
larized. This is done by cutting the spatial integrals off 
at the persistence length £ over which a polymer is stiff. 
This contains the stiffness caused by the excluded- volume 
effects. To be rigorous, we define the integral (|l5| ) on a 
lattice with spacing £. 

Replacing the expectation values by the Wick contrac- 
tions corresponding to the first diagram in Fig. ||, we 
obtain 



V M 4 



f dt[{l-t)t]-i j d 3 ye~^ J d 3 x f {^. 



xe ^(i-s) 



The variables x and y have been rescaled with respect to 
the original ones in order to extract the behavior of Ni 
in Li and L 2 . As a consequence, the lattices where x and 
y are defined have now spacings £/y/Li and ^j\fL 2 re- 
spectively. The x, y integrals may be explicitly computed 
by analytical methods in the physical limit L\^L 2 » £, 
in which the above spacings become small. This has a 
physical explanation. Indeed, if the polymer lengths are 



much larger than the persistence length, the effects due 
to the finite monomer size become negligible and can be 
ignored. 

Finally, it is possible to approximate the integral in x" 
with an integral over a continuous variable p and a cutoff 
in the ultraviolet region: / dV/l/lx" | 4 - 4tt 2 /~ dp/p 2 . 
After these approximations, we obtain 



(16) 



7. For the second diagram in Fig. || we have to calcu- 
late 

2 i. r +i °° Md Zl Mdz 2 t t 
Af 2 = ft 2 lim / ± e ziL 1+ z 2 L 2 

ZZlJc-ioo 2ttz 2m 
x J d 3 x 1 d 3 x 2 J d 3 x[d 3 x"d 3 x 2 

x^r(xi)| 2 |^2 2 (x2)| 2 (A 1 -^V^ 1 ) x , 

x(A 1 -^V^ 1 ) xr (A 2 |^ 2 | 2 ) x , y (17) 

The above amplitude has no ultraviolet divergence, so 
that no regularization is required. The Wick contractions 
pictured in the second Feynman diagrams of Fig. ^| lead 
to the integral 



i\/f3 ^1 ft 

N 2 = -A-sftVLz 1 ' 2 ^ 1 —^ / dt dt'C(t,t'), 
n Jo Jo 



(18) 



where C(t,t f ) is a function independent of L\ and L 2 : 



C(t,t') = [(l-t)t'(t-t')} 



/M-3/2 



d 3 xd 3 yd 3 ze 2 (' 1 - t ) 



with 

P^(x,,y,x) = [^ l/ z-(z+x)-(z + x) /i ^]/(|z| 3 |z + x| 3 ). 
As in the previous section, the variables x, y, z have been 
rescaled with respect to the original ones in order to ex- 
tract the behavior in L\. Again, if L\^L 2 » £ the an- 
alytical evaluation of C(t,t f ) becomes possible, leading 
to 



N 2 



vl: 



-l/2r-l 



(2tt)< 



^-M 3/2 4if, 



(19) 



where K is the constant \B (§, |) + |# (§, §) - 
B(l,\) + \B(%,\) = 19tt/384 « 0.154, and B(a,b) 
is the Beta function. For large L\ — > oo, this diagram 
gives a negligible contribution with respect to N±. 

The third diagram in Fig. || give the same as the sec- 
ond, except that L\ and L 2 are interchanged: A3 = 
N 2 \ Li ~l 2 . 

8. The fourth Feynman diagram in Fig. ^ has no ultra- 
violet divergence. As before, it can be exactly evaluated 
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apart from the lattice integrations. However, the be- 
havior of the related Feynman integral N4 can be easily 
estimated in the following limits: 

1. Li > 1; Li > P 2 , where 7V 4 oc L^ 1 , 

2. P 2 > 1; L 2 > Li where iV 4 oc P^ 1 , 

3. Pi,P 2 ^> 1, L2/L1 = finite, where ex P^ 2 . 
Moreover, if the lengths of the polymers are considerably 
larger than the persistence length, N4 can be computed 
in a closed form: 



Na 



128V M 



7T 7T 6 ' Z 



-1/2 



x f ds f dt(l-s)(l-t)(st) 1/2 
Jo Jo 

x [L 1 t(l-s)+L 2 (l-t)s}- 1/2 . 



(20) 



It is simple to check that this expression has exactly the 
above behaviors. 

9. Collecting all contributions we obtain the result for 
the second topological moment (m 2 ) = (Ni + N 2 + A3 + 
AM/Z, with iVi, 7V 2 , 7V 3 , 7V 4 , Z given by Eqs. @, ©, 
(|l9|), and (po|). In all formulas, we have assumed that the 
volume V of the system is much larger than the size of 
the volume occupied by a single polymer, i.e., V ^> L\. 

To discuss the physical content of the above expression 
for (m 2 ), we consider a number N p of polymers p\ . . - Pn p 
with lengths h . . .In p in an uniform solution. We intro- 
duce the polymer concentration p = M./V as the average 
mass density of the polymers per unit volume, where Ai 



is the total mass of the polymers M = 



N p 
k=l 



,1k/ a and 



m a is the mass of a single monomer of length a. Thus Ik /a 
is the number of monomers in the polymer pk. The poly- 
mer Pi is singled out as anyone of the polymers pk , say 
Pk, of length L\ = The remaining ones are replaced 



by a long effective polymer P 2 of length L 2 
From the above relations we may also write 



L 2 



aVp 



m a 



(21) 



In this way, the length of the effective molecule P 2 is 
expressed in terms of physical parameters. Keeping only 
the leading terms for V ^> 1, we find for the average 
square number of intersections (m 2 ) so i formed by P x with 
the other polymers the approximate result 



2 _ N! + N 2 



sol 



z ' 

which, in turn, has the approximate form 



(m 2 ) so i 



ap 

m a 



2KL 



1/2 



r 1 ^ 

27r l/2 M 2 7^3/2 



(22) 



(23) 



with K as defined after Eq. (19). This is the announced 
final result. Since the persistence length is of the same 
order of the monomer length a and M ~ a -1 , (m 2 ) is 
positive for large L\ as it should. 



10. In conclusion, we have found an exact field theo- 
retic formula for the second topological moment of two 
polymers. Only the final integrations over the spatial 
variables in the Feynman diagrams of Fig. || were done 
approximately. These were defined on a lattice related to 
the finite monomer size. Our Chern-Simons-based theory 
is free of the shortcomings of previous mean-field proce- 
dures. Our formula for (m 2 ) has been applied to the re- 
alistic case of long flexible polymers in a solution. When 
the polymer lengths become large, the Feynman integrals 
in Ni , . . . , N4 can be evaluated analytically. In this way 
we have been able to derive the result (|||) for the average 
square number of intersections formed by a polymer Pi 
with all the others. This calculation is exact in the long- 
polymer limit. The corrections to ( p3|) are suppressed by 
further inverse square roots of the polymer lengths. 

To leading order in Pi, our result ([23| ) agrees with 
that of H , but our exact subleading correction go beyond 
the approximation of p. Note that there is no direct 
comparison of our result with that of Q , since there the 
polymer P 2 was considered as a fixed obstacle causing a 
dependence on the choice of the configuration of P 2 . 

Finally, let us emphasize the absence of infrared di- 
vergences in the topological field theory ([I|) in the limit 
of vanishing masses m\ 1 m 2 =0. As a consequence, the 
second topological moment does not diverge in the limit 
of large Pi if (m 2 ) is calculated from (Q) for polymers 
passing through two fixed points xi,X2. This indicates 
a much stronger reduction of the configurational fluctu- 
ations by topological constraints than one might have 
anticipated. 
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